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Abstract 

To leading order in Ug , the leading and non- leading 1 / mb corrections to the excited 
B** meson coupling gB**Bn is calculated in the framework of QCD spectral moment 
sum rules in the full theory. Our prediction is in good agreement with the light — 
cone QCD sum rule result. 
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1 Introduction 



At present the physics of hadrons containing heavy quarks hes in the focus of experimen- 
tahsts and theoretists. From experimental point of view it is connected in the first hand, to 
the the observed exclusive B K*'y and inclusive decays B Xs'y 0, which clearly 
demonstrated the existence of FCNC and in turn, opens a window for the determination 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vtb and Vts- From theoretical 
point of view the reason is that the investigation of these decays allows one to check the pre- 
dictions of CKM at loop level. At the same time these decays form the basis of a very rich 
"laboratory framework" for theoretical researches. For example, in the infinite heavy quark 
mass limit (mg oo), QCD reveals symmetries that are not present in the finite mass 
theory: chiral SUl{3) ® SUr{3) symmetry, spin-flavour symmetry Qj. These symmetries 
allow the derivation of a simple relation among the formfactors of the different processes 
(see for review [^]) and build an effective heavy-light Lagrangian. This effective Lagrangian 
widely applied in the investigations of B and D meson physics. This approach allows us to 
include the positive parity {Qq) meson to the effective Lagrangian according to the value 
of the value of the angular momentum of the ligth degree of freedom (sf = ■^). The 
heavy quark effective theory p|, |[ predicts the existence of two multiplets, the first one 
containing 0+ and 1+ and the second 1"*" , 2+ states mesons. 

In ^ the strong coupling constant gB**Bn, where B** is the 0"^ p-wave in the framework 
of the classical and light cone sum rules, is investigated. Predictions of the both approaches 
differ one another by a factor ~ 1.5. This article is devoted to the study of the large 
(l/mfe) behavior of the gB**Bn coupling constants and estimate their values in the double 
momentum version of the QCD sum rules. Contrary to the double exponential version 
of sum rules, this version is advantageous in the analysis of the three point function, as it 
prevents the blow up of the QCD series, when mq is large. Note that the coupling constant 
gB*BTT^s calculated in the framework of the above mentioned method in [§]. 

2 Calculation of of the excited Qb^^b-k coupling 

We start our calculation by considering the following correlator: 

V = - jd'^x A e'^P'y-P-^\Q\T {JBA^).W)JB{y)} |0) , (1) 
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where the interpolating currents 



Jb** = nibuh , 
J5 = {rUu + nidjujsd and 

Jb = {md + mb)dj5b , (2) 

are described 0"*" excited B meson state in the sf — ^ doublet, 7r-meson and S-meson 
states. Acoording to QCD ideology this correlator is calculated theoretically in the deep- 
Eucledian region and in terms of physical states. In order perform this calculation we 
remind that the correlator obeys the double dispersion relation 

/ X 1 f°° ds ds' , „ 

V(j,,M = --J^^^—,l^^—^l^V(s,s') + ... (3) 

Since rub is much larger than the QCD scale A, it is natural to use the momentum ver- 
sion of sum rules. Acting on both sides of Eq. (3) by the double momentum operator, 
— f — over variables p^,p'^, and after differentiation setting — p'^ = 0, we get 

^(n,n') ^ _ 1 r ^ r Jf_ I^V{S, S') . (4) 

b b 

The perturbative part of the spectral function can be calculated from the Feynman diagrams 
by the help of the Gutkovsky rule, i.e., replacing propagators 

After carrying standard calculations for spectral density we get, 

1 Nc 
;lmV{s, s') = {mu + md)ml —- — — , (6) 



where = —q^ > is the pion momentum squared and 

X^{s + s' + Qy-4ss' . (7) 
Integration region defined by the following inequality: 

{s-ml){s'-ml)>Q'ml (8) 

The physical part of the sum rules is parametrized using the usual duality anzats: 
lowest resonances + QCD continuum, which usually modelled as bare loop, starting from 
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the thresholds Sq and Sg. Transferring the the continuum contribution to the theoretical 
part we obtain the following QCD sum rules 



93**371 



2(n+l) 



m 



2(n'+l) ^2 + Q2 4^2 ^n+l ^'n'+l 



*o (is /■''o ds' 



lmV{s,s'), (9) 



In deriving the left hand side of the sum rules we use the following definitions: 



9B**3tt 

{0\J,\n) 
{0\Jb\B) 
{Q\Jb**\B**) 



{Bip')n{q)\B**ip)) 



(10) 



In the irtb — >■ oo case it is convenient to work with the non-relativistic variables E and E' 
defined as 

s = {E + rribf and s' = {E' + rrikf , 
and introduce the new variables: 
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(11) 



x = E-E' and y = -{E + E') . 



(12) 



The integration region over new variables x and y is defined from Eq. (8). In further 
analysis we shall set n = n' = n^. By keeping the non-leading terms in the expansion to 
the leading order in a^: 



(m„ + rud) 



(x2 + Q 



dx I dv ~h / dx I 

2 

1 - (4n3 + 3)-^ + 2(^3 + l)(4n3 + 3)^ + 



2^1/2 



nib 



mt 



X 



^'""^^'Uml 8ml 



(13) 



Following 1^, H , we shall use in our analysis the lowest order expression in for the decay 
constants fs** and J'b from the 2-point momentum sum rule, 



27r2 m^VrnJ \ 2^' 'nib 2 (Ef)^J 



From the numerical analysis of fs including the correction one can obtain for the 
effective value of 

~ (1.3±0.1) GeV . (15) 

Performing similiar calculation using Eq. (14) and the results cited in for the effective 
value E^" of fs** we get, 

Ef" ~ (1.8±0.1) GeV . (16) 

Using the recently reported experimental results [|I^, [TTj] m^.. = (5.732 ± 0.002) MeV and 
rriB = 5.279 MeV we conclude that rriB** — itib = A ~ 0.5 GeV. From comparision of Eqs. 
(15) and (16) we have, 

Ef'^Ef + A. (17) 

Taking into consideration these values of E^" and m*^ and following |jl2|, we minimize 
the n-dependence of the results, by requiring that the leading term is n- independent. This 
leads us to the following condition: 



Am = 2n2 - 1 . 



From analysis of two point function it is found that ?22 — 4 — 5 (see [|13[ ) . Evaluating 
different integrals, the results can be written in the following good approximation, 

' En 



where: 



and 



LO 



Ncjmu + mrf)m^ 



Xn 



niB (^, 



B^3 



ax / ^ -, , 

J^y/x^+Q^ Jo 



dx 



rpB X 
Ec -2 



dy 



[x^ 



Q 



2^l/2 ■ 



(19) 



(20) 



(21) 



In the region 1 < < 4 this integral has a maximum at ~ 2.5 GeV^. Using this 
value of Q^, numerical integration of Xq for different values of Ej^ can be parametrized by 
following interpolating formula 



Jo = (0.175 ±0.001)Ec 



(22) 



Using 



E. 



B 



E': 



El 

2mh 



(see for example [8] 



(23) 
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we finally get 

Using the result g°° = (0.15 ± 0.03), that is given in ||^ and taking into consideration the 
correction for the physical B — meson (which is of the order of ±28% at the 6-mass), 
we get, 

gB**B^ = 24 ± 7 GeV , (25) 

where the error in the above result takes into account the effect of radiative corrections to 
Jb- As mentioned in Introduction, the value of the strong coupling constant gB**B-n that is 
calculated in |0| within the framework of the classical and light- cone sum rule techniques 
are , 

gB,.B^ = 13.3 ±4.8 GeV, 

gB**B. = 21 ±7 GeV, (26) 

respectively. From comparision of Eqs. (25) and (26) we see that our prediction on gB^Sw 
is in good agrrement with the light- cone QGD sum rule prediction on the same quantity. 
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